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It is shown that the quasilinearization method (QLM) sums the WKB series. The method ap- 
proaches solution of the Riccati equation (obtained by casting the Schrodinger equation in a nonlin- 
ear form) by approximating the nonlinear terms by a sequence of the linear ones, and is not based 
on the existence of a smallness parameter. Each p-th QLM iterate is expressible in a closed integral 
form. Its expansion in powers of h reproduces the structure of the WKB series generating an infinite 
number of the WKB terms. Coefficients of the first 2 P terms of the expansion are exact while coef- 
ficients of a similar number of the next terms are approximate. The quantization condition in any 
QLM iteration, including the first, leads to exact energies for many well known physical potentials 
such as the Coulomb, harmonic oscillator, Poschl- Teller, Hulthen, Hyleraas, Morse, Eckart etc. 

PACS numbers: 03.65.Ca, 03.65.Ge, 03.65.Sq 



I. INTRODUCTION 

The derivation of the WKB solution starts by casting 
the radial Schrodinger equation into nonlinear Riccati 
form and solving that equation by expansion in powers of 
h. It is interesting instead to solve this nonlinear equation 
with the help of the quasilinearization method (QLM) 
Q& H II II H and compare with the WKB results. 
The quasilinearization method and its iterations were 
constructed as a generalization of the Newton-Raphson 
method 0, Q to yield rapid quadratic and often mono- 
tonic convergence to the exact solution. It does not rely 
on the existence of any kind of smallness parameter. The 
initial comparison of QLM and WKB was performed in 
the work j^, EU , where it was shown that the first QLM 
iteration reproduces the structure of the WKB series gen- 
crating an infinite series of the WKB terms, but with 
different coefficients. Besides being a better approxima- 
tion than the usual WKB, the first QLM iteration is also 
expressible in a closed integral form. 

The goal of this work is to point out that the similar 
conclusions could be reached for higher QLM approxima- 
tions as well. Namely, we show that the p-th QLM iterate 
with p > 1 reproduce the WKB series exactly up to H 2 : 
when expanded in powers of h, it, besides providing the 
correct structure of the whole series, generates the coef- 
ficients of the first 2™ terms of the WKB series precisely 
and of a similar number of the next terms approximately. 
In addition we prove that the exact quantization condi- 
tion in any QLM iteration, including the first, leads to 
exact energies not only for the the Coulomb and har- 
monic oscillator potentials as it was shown in ref. 

HQS; 

but for many other well known physical potentials used in 
molecular and nuclear physics such as the Poschl- Teller, 
Hulthen, Hyleraas, Morse, Eckart etc. 

The paper is arranged as follows: in the second chapter 
we cast the radial Schrodinger equation into the nonlin- 
ear Riccati form and solve it by QLM iterations. Then in 
the third chapter we expand the first, second and third 
QLM iterates in the series in powers of H and and com- 
pare these expansions with the usual WKB series. In 



the next, forth chapter, we apply our method to equa- 
tions with the Coulomb, harmonic oscillator, Poschl- 
Teller, Hulthen, Hyleraas, Morse and Eckart potentials 
and prove that the QLM method provides exact energies 
for both ground and excited states already in the first it- 
eration. Our results, advantages of the quasilinearization 
method over WKB and its possible future applications 
are discussed in the final, fifth chapter. 



II. QLM APPROACH TO THE SOLUTION OF 
THE SCHRODINGER EQUATION 

The usual WKB substitution 

X (r)=Ce X p(\J r y(r')dr'y (1) 

converts the radial Schrodinger equation 

+ A 2 fc 2 X (r) = (2) 



dr 2 



to nonlinear Riccati form 



^l + (k 2 (z)+y 2 (z))=Q. 



(3) 



Here k 2 (z) = E - V - 1(1 + \)/z 2 , A 2 = 2m/h 2 and z = 
Ar. 

The proper bound state boundary condition for poten- 
tials falling off at z ~ zq ~ oo is y(z) — const at z > zq. 
This means that 

V'(zo) = 0, (4) 

so that Eq. © at z ~ zq reduces to k(zo) 2 + y 2 (zo)) = 
or y(zo)) = ±«fc(z ). We choose here to define the 
boundary condition with the plus sign, so that 

y(z Q ) = ik(z ). (5) 

The quasilinearization 0, 0, Q of this equation gives a 
set of recurrence differential equations 

= y l_ l{z) _ 2y p (z)y p _ l (z) - k 2 (z). (6) 
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with the boundary condition deduced from Eq. ((SJ): 

y P (z a ) = ik(z ). (7) 

The analytic solution |9j of these equations expresses 
the p-th iterate y P {z) in terms of the previous iterate: 



Vp{z) = fp-i( z ) 



where 



f P -i( z ) 



df p -i(s) 

as — — exp[— 2 

as 



y p -i(t)dt], 
(8) 

(9) 



1 dfp^js) 
2j/ p _i(s) ds 



Indeed, differentiation of both parts of Eq. © leads im- 
mediately to Eq. jnj which proves that y p (z) is a solution 
of this equation. The boundary condition JJJ is obviously 
satisfied automatically. 

The second term in Eq. (JSJ could be written as 



ds 



dfp-l(s) 



2y P -i(«) 



2y p _i(s) exp[-2 / y p -i(t) dt] 



(10) 



The second expression in the round brackets in the in- 
tegrand is the derivative of the exponential there. The 
integration by parts of this integral therefore gives: 



exp[-2 f y p -i(t)] - f 



ds- 



dfp-l(s) 
ds 



ds \ 2y p -i(s) 



exp[-2 / y p ^(t))dt], 



(11) 



Since the lower limit in the first term of this expression term of Eq. JHJ by parts results in: 
vanishes in view of Eq. 10} , the integration of the second 



y P i z ) = f P -i( z ) + 



df p -i{z) 



2y p _i(z) dz 



pz l I rf /p-l( s ) \ rz 



(12) 



The successive integrations by parts of Eq. (|12|l lead 
to the series 



(13) 



compare with the WKB series as will be done in the next 
section. 

For the zeroth iterate yo(z) it seems natural to choose 
the zero WKB approximation that is to set 



with C^n\z) given by recursive relation 

C%\z) = —^±(-C%Uz)) 
2y p -i(z) dz 



y (z) = ik(z), 



(17) 



and 



Since 



4 0) « = / P -iW- 



— = — =A~ 1 
dz dr ' 



'2m ' 



(14) 



(15) 



(16) 



Eq. (|13fl represents the expansion of the p-th QLM iter- 
ate in powers of <?, that is in powers of h, which one can 



which in addition satisfies boundary condition JSJ. 
However, one has to be aware that this choice has 
unphysical turning point singularities. According to the 
existence theorem for linear differential equations if 
y p ^i(z) in Eq. © is a discontinuous function of z in a 
certain interval, then y p (z) or its derivatives may also be 
discontinuous functions in this interval, so consequently 
the turning point singularities of yo(z) may propagate 
to the next iterates. 

Substitution of the initial QLM approximation l|17l) 
into Eq. (JSJ) gives especially simple expression [9j for the 



3 



first QLM iterate 



yi{z) = ik{z) — i / dsk'(s) exp[— 2i / k(t)dt], 



(18) 

Thus the first QLM iterate is expressible in a closed 
integral form. However, it takes into account, though 
approximately, infinite number of the WKB terms corre- 
sponding to higher powers of h as it will be shown in the 
next section. In view of this it is a better approximation 
than the usual WKB. 



III. COMPARISON OF EXPANSIONS OF QLM 
ITERATES AND WKB SERIES 

To obtain the WKB series one has to expand solution 
y of the Riccati equation (JHJ in powers of fi.This is easy 
to do by using Eq. lltjfl and looking for y in the form of 
series in g: 



E 



iy 

- 1 'II, 



(19) 



Substitution into Q and equation of terms by the iden- 
tical powers of g gives 



dr 



^ Yk Y m —k- 



(20) 



This reduces to the recurrence relation 



1 / 

Ym = ~2% y™- 1 + E 



Y k Y 



m—k 



(21) 



The derivatives in this and subsequent expressions are 
in variable r. The zero WKB approximation Y Q is given 
by Yq — i k. The subsequent terms Y n of the expansion 
could be obtained from this recurrence relation by use of 
Mathematica 

We present here the WKB expansion (|19|) up to g 7 
inclusively: 



y = z fc - ^ + 4_ 5 2 (3 fc' 2 - 2 fc k") + -4= (6 k' 3 - 6 k k' k" + k 2 + —!— g 4 (-297 fc' 4 + 396 k k' 2 k" 
2 k 8 fc 3 y ' 8 k 5 v ; 128 k 7 y 

5 

-52 k 2 k" 2 - 80 k 2 k' k® + 8 k 3 k^) - J— (306 fc' 5 - 510 fc fc' 3 fc" + 111 fc 2 fc' 2 k ( ^ _ 3 fc 2 fc' (-48 fc" 2 + 5 fc fc< 4 >) 

32 fc y 

+fc 3 (-24 fc" fc( 3 > + fc fc( 5 >)) + 1Q2 ^ fcll 9 e (50139 fc' 6 - 100278 fc fc' 4 fc" + 22704 fc 2 fc' 3 fc (3) + 12 fc 2 fc' 2 (3679 fc" 2 
-290 fc fc (4) ) + 16 fc 3 fc' (-694 fc" fc (3) + 21 fc fc (5) ) - 8 fc 3 (301 fc" 3 - 80 fc fc" fc (4) + fc (-49 fc (3)2 + 2 fc fc (6) ))) 

7 

— ^—r (38286 fc' 7 - 89334 fc fc' 5 fc" + 20721 fc 2 fc' 4 fc (3) + fc' 3 (53724 fc 2 fc" 2 - 3405 fc 3 fc (4) ) + 3 fc 3 fc' 2 (-5426 fc" fc (3) 
128 k Li 

+ 129fcfc (5) ) + 2fc 3 fc'(-3528fc" 3 + 735fcfc"fc (4) +2fc(225fc (3)2 - 7 fc fc (6) )) + fc 4 (1176 fc" 2 fc (3) - 62fcfc"fc (5) 

+fc(-90fc (3) fc (4) +fcfc (7) ))) 

(22) 



To compare expansion of the first QLM iterate y\ in pow- 
ers of h with the WKB expansion l|22|l we have to use, as 
we have already mentioned in the previous section, Eqs. 



(H2J and (H1J) together with Eq. JT7J)- The result up to 
power g 7 inclusively is again obtained with the help of 
Mathematica [lj: 
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yi = ik ~ + ^ (k' 2 - kk") + ^{3k' 3 - Akk' k" + k 2 k<-V ) + (_15 fc' 4 + 25 fc fc' 2 fc" - 7 fc 2 fc' fc' 3 ) 

+fc 2 (-4 fc" 2 + fc fcW)) + -^—(-105 fc' 5 + 210 fc fc' 3 fc" - 60 fc 2 fc' 2 fc^ + fc 2 fc' (-70 fc" 2 + 11 fc k^) 

32 fc y 

+fc 3 (15 fc" fc (3) - fc fc (5) )) - (-945 fc' 6 + 2205 fc fc' 4 fc" - 630 fc 2 fc' 3 fc (3) + 14 fc 2 fc' 2 (-80 fc" 2 + 9 fc fc (4) ) 

64 fc 11 

7 

+2 fc 3 fc' (175 fc" fc (3) - 8 fc fc (5) ) + fc 3 (70 fc" 3 - 26 fc fc" fc (4) + fc (-15 fc (3)2 + fc fc (6) ))) + g (10395 fc' 7 

128 K 

-27720 fcfc' 5 fc" + 7875fc 2 fc' 4 fc (3) - 126fc 2 fc' 3 (-150 fc" 2 + 13fcfc (4) ) + 14fc 3 fc' 2 (-495fc"fc (3) + 17fcfc (5) ) 
+fc 3 fc' (-2800 fc" 3 + 784 fc fc" fc (4) + fc (455 fc (3)2 - 22 fc fc (6) )) + fc 4 (560 fc" 2 fc (3) - 42 fc fc" fc (5) 

+fc(-56fc (3) fc (4) +fcfc (7) ))) 

I 



The comparison of expansion of the first QLM iterate in 
powers of h and WKB series was originally performed in 
the works 0, 0] . There it was shown that the expansion 
reproduces exactly the first two terms and also gives cor- 
rectly the structure of the WKB series up to the power 
g 3 considered in these works, generating series with iden- 
tical terms, but with different coefficients. Comparison 
of Eqs. (|22|l and \23\ ) of the present work shows that this 



conclusion is true also if higher powers of g are taken into 
account. 

The computation of the expansion of the second QLM 
iterate 1/2 in powers of H is done by reexpanding the term 
2^- in Eq. (|14|l in the series in powers of g and keeping 

the powers up to g 7 inclusively in this expression as well 
as in the sum in Eq. (|13|) . The result is given by: 



y 2 = ifc-^# + ^(3fc' 2 -2fcfc") + -^(6fc' 3 -6fcfc' fc" + fc 2 fc^) + ^(-74fc' 4 + 99fcfc' 2 fc"-20fc 2 fc'fc^ 
2 fc 8 k 6 8 fc° 32 fc' 

+fc 2 (-13 fc" 2 + 2 fc fc( 4 >)) + -ija (- 607 k ' 5 + 1017 fc fc' 3 fc" - 222 fc 2 fc' 2 fc( 3 ^ + 6 fc 2 fc' (-48 fc" 2 + 5 fc fc (4) ) 

64 fc y 

-2 fc 3 (-24 fc" fc (3) + fc fc (5) )) - if (-6186 fc' 6 + 12446 fc fc' 4 fc" - 2832 fc 2 fc' 3 fc (3) + fc 2 fc' 2 (-5503 fc" 2 + 435 fc fc (4) ) 

7 

+2 fc 3 fc' (694 fc" fc (3) - 21 fc fc (5) ) + fc 3 (301 fc" 3 - 80 fc fc" fc (4) + fc (-49 fc (3)2 + 2 fc fc (6) ))) + 9 1Q (75256 fc' 7 

256 k Li 

-176659 fc fc' 5 fc" + 41224 fc 2 fc' 4 fc (3) + 4 fc 2 fc' 3 (26687 fc" 2 - 1700 fc fc (4) ) + 2 fc 3 fc' 2 (-16243 fc" fc (3) + 387 fc fc (5) ) 
+fc 3 fc' (-14071 fc" 3 + 2940fcfc"fc (4) +8fc(225fc (3)2 - 7 fc fc (6) )) + 2 fc 4 (1176 fc" 2 fc (3) - 62fcfc"fc (5) 

+fc(-90fc (3) fc (4) +fcfc (7) ))) 

(24) 



The expansion of yi reproduces exactly already the first 
four terms of the WKB series. It also gives the proper 
structure of the other terms of the WKB series, generat- 
ing series with identical terms which have approximately 
correct coefficients. 

The expansion of j/3 is obtained in the similar fashion. 
It reproduces exactly the first eight terms of the WKB 
series that is all the terms up to the power g 7 inclusively 
listed in Eq. JUJ. 

Summing up, we have proved that the expansion of 



the first, second and third QLM iterates reproduces ex- 
actly two, four and eight WKB terms respectively. Since 
the zero QLM iterate yo was chosen to be equal to the 
zero WKB approximation i fc, one can state that the p-th 
QLM iterate contains 2 P exact terms. In addition expan- 
sion of each QLM iterate has the correct structure whose 
terms are identical to of the WKB series with approxi- 
mate coefficients. 

The 2 P law is, of course, not accidental. The QLM 
iterates are quadratically convergent 0, 0, E| , that is the 
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norm of the difference of the exact solution and the p-th 
QLM iterate \\y — y p \\ is proportional to the square of 
the norm of the differences of the exact solution and the 
(p-l)-th QLM iterate: 



\\y-v P \ 



\y-y P -i\ 



(25) 



Here norm || g || of function g(x) is a maximum of the 
function g(x) on the whole interval of values oi x. 

Since yo contains just one correct WKB term of power 
g° and thus \\y— yo\\ is proportional to g, one has to expect 
that || y— 2/1 1| ~ g 2 and thus y\ contains two correct WKB 
terms of powers g° and g . The difference \\y — J/2II ~ 
1 1 2/ — J/x I j ~ g 4 so that 2/2 contains four correct WKB 
terms of powers g°, g , g 2 and g 3 . Finally, the difference 
\\y — 2/3 1| should be proportional to g s , and therefore y^ 
has to contain eight correct terms with powers between 
g° and g 7 inclusively. This explains the 2 P law. 



IV. QLM AND WKB ENERGY CALCULATIONS 

The exact quantum mechanical quantization condition 
for the energy 0,0,0] has the form: 



y{z)dz = i 2nn. 



Here y(z) is the logarithmic derivative of the wave func- 
tion, given by Eq. ©, z — gr, n = 0, 1, 2, ... counts the 
number of poles of y(z) and is the bound state number 
and the integration is along a path C in the complex 
plane encircling the segment of the %iz axis between the 
turning points. 

The p-th QLM iterate y p (z), as we have seen, con- 
tains in addition to 2 P exact WKB terms of powers 
also an infinite number of structurally 



■ ff 



2 P -1 



correct WKB terms of higher powers of g with approx- 
imate coefficients. One can expect therefore that the 
quantization condition (|26|l with y(z) approximated by 
any QLM iterate y p (z) including the first 



J p = j> y p (z)dz = i 27rn, p = 1, 2, ... (27) 

gives more accurate energy than the WKB quantization 
condition 



k(z)dz = 2tt(ti + -). 
c 2 



(28) 



which is obtained by substituting into exact quantization 
condition (gSJ the WKB expansion 

J23 up to the first 

dk(z) 



2k{z) 



and ne- 



power of g ~ H, that is y(z) = ik(z) 
glecting all higher powers of g in the expansion. Indeed, 
we prove in this section that Eq. (|27|) leads to exact 
energies not only for the the Coulomb and harmonic os- 
cillator potentials as it was shown earlier in ref. 0, ^| , 



but for many other well known physical potentials used in 
molecular and nuclear physics such as the Poschl- Teller, 
Hulthen, Hyleraas, Morse, Eckart etc. The WKB quan- 
tization condition 128|) yields the exact energy only for 
the first two potentials, but not for the rest of them. 
Let us now consider the above mentioned examples: 



A. Harmonic oscillator V(x) 



-00 < x < 00. 



1 so 



From now on we will work in the units h — m 
that from z — Xx follows x = and V(z) = 

In view of the boundary condition (JJJ y p (z) at the in- 

— where we 



finity should behave like i^j E — ^- ~ — | 

omitted terms of order \ and higher. Here we took into 
account that for bound states the logarithmic derivative 
at the infinity should be negative. More accurately the 
pole structure of y p (z) at z ~ 00 could be found by look- 
ing for the solution there in the form y p (z) ~ — -| + 
Substituting into the quasilinearized equation JHJ 

^ = y^z) - 2yp{z)yp„ l {z) - (E - ^ )) (29) 



and again neglecting terms of order and higher which 
(26) does not contribute to the integral yields a p — E — |, so 

E-i 



that the pole term in y p (z), p = 1,2, ... is given by 

The integration in Eq. I|27|) is counterclockwise along 
a path C in the complex plane encircling the segments 
of the $tz axis between the two turning points —2\[E 
and 2-y/E . Since the only singularity outside contour C 
in the complex plane lays at infinity, the integral 127|) 
could be done by distorting the contour to enclose the 
pole at x — 00 . The evaluation of the integral yields 
i2n (E — I ) = i2-Kn or E = n+ h, which is the exact 
equation for the energy levels. 



B. Spherical harmonic oscillator 



V(r) 



< r < 00. 



This case was already considered in the work 10]. We 
present here somewhat different derivation. Eq. © now 
has the form 



dy P {z) 
dz 



y p -i(z) - 2y p (z)yp- 1 (z) -(£- — - 



z 2 1(1 + 1) 



)) 



(30) 

where z = \[2r . The pole structure at z ~ 00 as before 
will be given by Eq. H3U|) . In order to find the pole 
at z ~ we look for y p near zero in the form -£ + .... 
where terms proportional to -4, etc. are not explicitly 
displayed since they do not contribute to the integral. 
The substitution into Eq. l|3Ufl yields 



1(1 + 1) 



(31) 
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When p becames large we expect that QLM iterates con- 
verge to each other and to the exact solution, so for large 
p one expects a p -\ = a p . The solution of Eq. (pTTJ) satis- 
fying this condition will be a p = (I + 1), so the pole term 
at z ~ is This, of course, agree also with a z l+1 
behavior of the radial wave function at zero, which leads 
to the behavior of its logarithmic derivative. The 
contour C encloses two positive turning points Z\ and z-i 
defined by the equation 



E 

4 



z 2 1(1 + 1) 



= 



(32) 



and the section 3?z axis between them. Indeed, the phys- 
ical motion takes place only for positive real z. The 
integrand, as we saw, has the poles at z = and 
z — oo. However, in addition 0, 0], when C is dis- 
torted to enclose these poles there will be also contri- 
bution to the integral J p from the cut on the nega- 
tive -Kz axis between the unphysical turning points — z-i 
and —Z\. This contribution due to the symmetry of 
the potential with respect to z — will be — J p since 
the path of integration around the cut on negative real 
axis will be clockwise. The contribution of the poles is 
2tti multiplied by — (I + 1) and E — i respectively, so 
one obtains J p = —J p + 2m[E — -| — (I + 1)]. Thus 
J p = m(E - i - (I + 1)) = 2inn or 



E = 2n + I 



(33) 



which is the exact equation for the energy levels of the 
three-dimensional oscillator. 



previous case there is no cut along the negative part 
of the real axis and when C is distorted to enclose 
the poles at zero and the infinity only contribution of 
these poles should be taken into account. One obtains 
2ixi(—^ (1 + I)) = i2im or E = ft*, isa which 

\/2\E\ v 2(n+i+l) 2 

are exact energy levels in the Coulomb potential. 



D. Cotangent potential 

V(x) = V cot 2 (?f), Vo>0, 0<x<a. 



Let us introduce a new variable 

• 2 / 7ra; \ 
z = sin ( — ) 

a 



(34) 



so that x 



zarcsiriy 



dx 



2tt 



dz and k(z) 



2(E + Vq) - ^. The QLM equation P will now have 
a form 



a dz 



V 2 p -i{z) - 2y p (z)y p -i(z) - k 2 (z). 

(35) 



One of the singularities of k 2 (z) is at z = 0. Near this 
point the equation has a form 

2?r r dy p (z) 2 2V . . 

— Vz ^ = Vp-iiz) - 2y p (z)y p - 1 (z) + —. (36) 

We look for solution of this equation in a form y p (z) = 
. Then we obtain for a p the following recurrence rela- 
tion: 



C. Coulomb potential V(r) 



< r < oo. 



This case was considered before in the work Q • In our 
derivation we set, as in the previous section z — v2r, 
so k 2 (z) = E + ^ - -^±±1 . At large z, in view of 
the boundary condition J^) y p (z) at the infinity should 

z 



behave like 



E 



Z^2 



= . where we 



took into account that for bound states the logarithmic 
derivative at the infinity should be negative and omitted 
terms of order \ and higher. The residue of the pole 
at the infinity thus is , z . The residue of the pole at 



z ~ is computed as in the previous section and yields 
the same value I + 1. 

The contour C encloses two turning points defined by 



and the section 



the equation -\E\ - ^ - ^11 
3tz axis between them. This equation, unlike equation 
(|32|l for the turning points in the previous section which 
had two positive and two negative real roots, has only 

z / z 2 i{i+T) 

2|^y yjw ~wr 

Therefore unlike the 



two positive real roots z\ = 
and zi - 




a p {2a p - x --) = a 2 p _ Y + (-) 2 A(A - 1). 



(37) 



Here we set 2Vq 



2V a 2 



(f ) 2 A(A - 1) where A 



The solution of this algebraic equation 
which fullfils the demand that at large p a p = a p _i is 
a p = —A. The y P (z) near zero thus has a form y p (z) — 
At z ~ oo Eq. EH reduces to ^y/^^^l = 

vl-ii z ) - 2yp( z )y P -i{ z ) - 2( E + V ). Looking for so- 
lution in a form y p (z) = c p where c p is some con- 
stant, one obtains the recurrent relation for c p , namely 
c 2 _ 1 - 2c p c p -i — 2(E + Vq) = 0. The solution of this al- 
gebraic equation which fullfils the demand that at large 
pc p = c p _i is c p = tJ-2(E + Vq). 

The quantization condition (|27(l in variable z given by 
Eq. has a form 



d,i 



y P (z) 



C y/z{l - Z) 



dz — i 2im, p — 1,2, 



(38) 



where integration is counterclockwise along a path C in 
the complex plane encircling the cut along the 5Rz axis 
between the zero and 



Vb 



2(E+V ) 



Since y p (z) equals to 
and y/—2(E + Vq) at z ~ and z ~ oo respectively, 
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the integrand in Eq. l(3*5jl has poles with residues ^ and 
\J2{E + Vq) there. The deformation of the contour to in- 
clude these poles and computation of their contributions 
yields 



J p = 2m{ A + y/2(E + V a )) = 2mn, 



(39) 



or, upon substitution of value of A and inserting h and 
?7i from dimensional considerations, 



E = -Vn 



^h" I , 1 , / 2mV , In 
2^aZ[ n + 2 + V tt 2 /1 2 + V' 

This is the exact equation for the energy levels in the 
cotangent potential The correspondent WKB ex- 

pression is different and given by 



E 



- V + K n + 2 + V 7T 2 /i 2 J 



E. Poschl- Teller potential hole 



Vl 

sin 2 (-^) 



+ 



V 2 



Vi > 0, V 2 > 0, < x < §. 



This potential is a generalization of the cotangent po- 
tential and reduces to it in case Vi = Vi = Indeed, in 
this case V 



Vn 



M«>i 2 (D + 1] where y = 2x 
The computation therefore 



sin 2 ( 

changes between and a 
proceeds as in the previous section with a difference that 
in addition to the poles at zero and the infinity there is 
also a pole at z = 1. Indeed, using again variable z de- 
fined by Eq. (|34|) one writes the Eq. (JSJ in a form i|35|) 
with k 2 (z) now given by 2(E - - j^). 

The poles at z = 0, z = 1 and z — oo are computed 
as in the previous section and have respectively the form 

f2E 



~7rz~ i t^tt and 



Here A*. = 1 
- 1); k 



i + W so 



2 1 V 4 

1,2. The integration in 



az 7 a(z — 1) 

that 2V k = (f ) 2 A fc (A fe 
the quantization condition (|27l is counterclockwise along 
a path C in the complex plane encircling the cut along 
the 3?z axis between the two turning points z\ , Z2 given 
by 21,2 = 5i[0E + Vi - V 2 ) ± V(^ + Vi-V 2 ) 2 - Wi] . 
The deformation of the contour to include the poles and 
computation of their contributions yields 



J p = 2tu(--Ai - -A 2 + V2E) = 2nin, (40) 

7T 7T 



or, upon substitution of values of Afe and inserting h and 
m from dimensional considerations, 

^h 2 



E 



1 



2m Via 2 
7r 2 R 2 



2roV2a 2 
ir 2 ri 2 



1)2 



2ma 2 Ln 

This is the exact equation for the energy levels in the 
Poschl- Teller potential hole 0|- The correspondent 
WKB expression is different and given by 



E : 



M2r 



1 



7T 2 /i 2 



2mV 2 a 2 l2 
7r 2 fi 2 J • 



F. Modified Poschl- Teller potential 



V ( x ) = - C os$(^) > Vo > 0, -oo < X < oo. 



Setting z = cosh 2 (-) so that dx — „ , a 

G n a ' 2z(z—l) 

tains that the QLM equation ijfJJl now has a form 

dy p (z) 



dz one ob- 



-V(z^z 
a 



(41) 

where k 2 (z) is given by 2(-\E\ + At z ~ this 

equation has a form 

2n dy p (z) 



dz 



y P -i( z ) - 2 y P (z)y P -i(z) 



A(A-l) 
Z (42) 



where we use the definition ofA = j; + y j + 2Voa 2 , so 

that 2Vq — ). The solution near zero could be 

looked for in analogy with previous cases in the form 



y P (z) 



Substitution into l|42|l gives 



Op-l 



2a p a p _i+A(A— 1) whose solution, satisfying the condition 
that at large p will be a p — a p -\ is a p = A — 1. Thus near 
zero y P {z) ~ r} . The solution at the infinity y p (z) ~ 



■\/2|_B| is found in the same way as in the previous two 
sections. The integration in the quantization condition 



Jp = a 



y P (z) 



C y/z(z - 1) 



dz = i 2nn, p = 1, 2, . 



(43) 



is counterclockwise along a path C in the complex plane 
encircling the cut along the dtz axis between z — = 

X 2a 2 \E\ anc ^ z = oo. The deformation of the contour and 
computation of the integral upon substitution of values 
of A and insertion of H and m from dimensional consid- 
erations yields 



E = - 



2ma 2 



-(n+-) + 



n 2 



(44) 



This is the exact equation for the energy levels in the 
Poschl-Teller potential hole (lHUri The correspondent 
WKB expression is different and given by 0] 

(n+ i ) + v sz ]2 ^ 



E 



2ma 2 I 



G. Hylleraas potential 

H*) = V > 0, 0<r <& o. 

This potential, used in the molecular and nuclear 
physics, was introduced in work 18] . The computation 
proceeds as in previous section with the difference that 
now the physical motion takes place only for positive real 
r and the wave function equals zero at the origin. The 
latter in view of the simmetry of the Hylleraas poten- 
tial toward exchange r to — r means that the solution of 
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the present problem is equivalent to considering uneven 
states in the potential of the previous section. Thus in 
Eq. (@3J n should be changed to (2n — l),n = 1,2,.. 
.Therefore the energy levels in the Hylleraas potential 
are given by 



ti 2 1 
2ma 21 y 2' 



2rnVoa 2 

h 2 



(45) 



This coincides with the exact expression for the 
energy levels. The correspondent WKB expression is 
rather different 



E 



=1,2, 



27ri(v^e + \J A — e — y — B) which yields to the following 
equation for the WKB energy levels: 



V2 



— (yfe + V^A - y/B) = (n + -) 



(48) 



This equation coincides with given in work |l9j . 

To compute energy levels in the quasilinear approxi- 
mation we have to use the QLM equation © which after 
switching to variable t has a form 



t dy 



B 



(I-*) 2 



(49) 



H. Eckart potential 



V(x) 



A- 



B- 



(l + e 



B > \A\, — oo < x < oo. 



This potential, introduced in work [20j, like the Hyller- 
aas potential is also widely used in molecular and nu- 
clear physics. The potential V(x) — ► —A as x — > — oo 
and V(x) — > as i — > oo. The minimum value 

V^i„ of this potential is V m i n = — 4^ ; the dis- 
crete spectrum of energy lies respectively in the interval 

(-i*+>pl,min(0,-A)). 

Let us start from calculating the energy levels e = — E 
in the WKB approximation. The WKB quantization 
condition (|28ll now has the form 



/2m / 
—f 



A- 



B- 



c 



(1 



-dx 



= 2n(n + -). (46) 

Introducing a new variable t = — e - ^,— 00 < t < re- 
duces the WKB quantization condition (|46ll to 




-e- A- 



B- 



"1-t (1-t) 2 t 



= 2n(n+-). (47) 

The contour C encloses two turning points t% t 2 = 
1 



2(e-A) 

equation e 



[(A + .B - 2e) ± v/(A + B) 2 - 4eB] defined by the 



A i-t 



B- 



= and the section 3?t 



of the axis between them. When C is distorted it en- 
closes the poles at zero, at one and the infinity so the 
contributions of these poles with residues y/^e, \/—B 
and — VA — e, respectively, should be taken into account. 
Since the poles at t = and t — 1 are enencircled clock- 
wise that is in opposite direction compare with the pole 
at t = 00 which is enencircled counterclockwise, their 
contribution enters with opposite signs and equals to 



For convenience of further computations we set here jfi 
equal to unity. The quantization condition (|27l) in vari- 
able t is given by 

J p = -a(f ^^dt = i 27171, ;p=l,2,... (50) 
J c t 

At the singular point t ~ of the integrand of Eq. I|5U|) 
Eq. ({4*9")l reduces to 



t dy p 2 



(51) 



whose solution is y p = c p , where c p is a constant sat- 
isfying the algebraic equations c 2 — 2c p c p _i + e. Since 
at large p we expect y p -\ — ► y p — > y, where y is an 
exact solution at t = 0, it means , in view of c p being 
a constant, that should be c p — c p ~\ — c, that is we 
are looking for a fixed point solution of Eq. 151(1 which is 
c p = The positive sign before the root is chosen since 
the first term in expansion of y P (t) in the WKB terms is 

ik(t). Thus y p (t) ~ ik(t) = ^Je + A-±- t + Bj-^—i that 

is y p (0) ~ +Ve- 

Another singular point of the integrand of Eq. 1(501) 

lays at t ~ 00 , since change of variable u = | , t — > 00 

when u 



0, converts 



eft into 



y P i v ) 



dv which has 



pole at v — 0. At t ~ 00 Eq. ((49(1 reduces to 



^=^-2j /p j /p _ 1 + e-A 



(52) 



The solution of this equation one can look in the form 
Up = d p where d p is a constant which satisfies the alge- 
braic equation d 2 — 2d p d p _ \ + e — A = 0. The fixed point 
d p = d of this equation is d = ye — A so that the residue 
of the integrand of Eq. ((50(1 at t ~ 00 equals ve — A. 
Eq. 1(49(1 near its singular point f = 1 has a form 



1 dy p 2 



A- 



1 -t 



B 



(i-ty 



(53) 
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Looking for solution of this equation in the form y p = j^r 
we obtain for constants b„ the recurrence relations -b„ = 
bp — 26 p 6p_i + B whose solution at the fixed point b p = b 
of this equation is given by b = — Vl + 4a 2 B) 

and y p ~ at t ~ 1. Again, the positive sign before 
the root is chosen from the same consideration as in the 
previous paragraphs namely since in zero WKB approx- 
imation the residue of y P (t) at t ~ 1 should be of order 
+VB. 

Summing up, the integrand of Eq. I|5(J|I has the three 
residues y/e, \Je — A and — ^(1 — \Jl + Aa 2 B) at t = 0, 
t = oo and t = 1 respectively. After the reinstatement of 
the factor ^ and taking into account that poles t = 
and = 1 are encircled in the opposite direction compare 
with the pole at the infinity (which is encircled counter- 
clockwise) the expression (|5U|1 gives therefore 



2ma 2 



1 



A) --(I 



ima 2 B 



(54) 



This expression coincides with the exact expression for 
the Eckart potential given in ref. |19| and is different 
from the WKB quantization condition (|48|l for the energy 
obtained earlier in this section. 



2ma 2 
If 



(-Ve + Ve + A) 



8m6 
If 



(57) 

which is different from the WKB expression i|56|l and 
coincides with the exact expression for the energy levels, 
calculated in work [23j |. 



J. Three dimensional S-wave Hulthen potential 



V(x) 



-X- 



A > 0, < r < oo. 



This potential is used in the nuclear physics and is a 
special case of the Eckart potential with 6 = 0. The 
Eqs. (jSEl) and for the WKB and QLM energy levels 
respectively degenerate to 



2ma 2 



(— v/e + Ve + A) = n + 



(58) 



and 



I. Three dimensional S-wave Eckart potential 



V(x) 



-X- 



+ b- 



^T 2 



A, b > 0, < r < oo. 



This potential is considered in works |21ll22ll23| where 
exact and WKB expressions for the energy levels were 
obtained. Introduction of a new variable t — e _ ",0 < 
t < 1 reduces the WKB quantization condition l|28l) to 



2ma 2 



(59) 



so that the WKB and QLM energy eigenvalues are ex- 
plicitly given by somewhat different expressions 



h 2 {{n- 



1)2 



2ma X \ 2 



2m 4a 2 (n+i) 2 



(60) 




-e+ A 



dt 



1 - t 



(l-t) 2 t 



2n(n+-). (55) 



The integrand has poles at t = 0, t 
oo with the residues V^e, y— b and v 



1 and t 



e — A, respec- 
tively. The calculation of the integral therefore gives |23| 
2ni(— \f^e — ^/— b + \f—k — A so that the WKB energy 
levels could be computed from expression 



2ma 2 
If 



f-Ve + Ve + A - Vb) = n + 



(56) 



On the other side, computation of the energy levels in the 
quasilinear approximation, following the steps, outlined 
in previous sections leads to the expression 



and 



h 2 ((n 



1\2 2ma A\2 

l > TP / 



2m 



Aa 2 {n + 1) 



(61) 



Here n = 0, 1, 2, ... The last expression coincides with 
the exact exp ression for the energy levels calculated in 
works [H[li|. 



K. One dimensional Morse potential 

V(x) = Ae~ 2 5 - Be"S , A,B > 0, -oo < x < oo. 

This potential, introduced in work [24j is heavily used 
in molecular physics computations to describe interac- 
tions between two molecules and for description of vi- 
brations of two-atomic molecules. Introduction of a new 
variable i = e",0<i<oo changes the WKB quantiza- 
tion condition l|28|) to 
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V2r 



A 



B dt 

7 T 

L 



2n(n+-) 



(62) 



The integrand has poles at t — oo and at t = with the 



residues 



and 



2V-A 

the integral therefore gives |2 



respectively. The calculation of 



V 2ma 2 , r- B . . 1 . 
— — _) = n+- 

which yields the well known expression 0, [2^, 



(63) 



The solution of this equation near t = has a form 
j/ p = ^f- + d p where a p and 6 P are constants satisfying the 
algebraic equations 



2bpb p -i 



A = 0, 



2bpd p — 2(6 p dp_i + bp-idp) - B. 



(69) 
(70) 



The fixed points b p = and d p — d p _i = d of this 

equation are b — \f~A and d = From all this 

the residues of the integrand at t = oo and t = are 
Y^e and ^ — ^= respectively. Taking into account that 
the pole at t = is encircled in the opposite direction 
compare with the pole at the infinity (which is encircled 
counterclockwise) the expression (|65H gives therefore 



(64) 



for the energy levels. 

To compute energy levels in the quasilinear approxi- 
mation one has to use the QLM equation JJjJ which in 
variable t has a form 



t dy p 
a dt 



yl - 2y p y p ~i + e + 



.4 



B 

T 



(65) 



We set here and further equal to unity. The quanti- 
zation condition l|27|l in variable t is given by 



yp(t) 



dt = i 2nn, p = 1, 2, . 



(66) 



a(7i-(-L-_^) 
2a 2VA 



(71) 



This expression after reinstating the factor ^™ coincides 
with H63|) which means that the QLM expression for the 
energy levels in the Morse potential coincides with the ex- 
pression H64fl for the WKB energy levels obtained earlier 
in this section. It coincides also with the exact expres- 
sion for the energy levels, see, for example, ref. |25j. The 
coincidence of the exact and the WKB eigenvalues for 
the Morse potential is consequence of the fact that the 
introduction of a variable u = 2a ^^ mA e - f ; q < y < oo 
reduces the Schrodinger equation with this potential to 
the radial Coulomb Schrodinger equation |25| which, as 
it is well known, yields exact energy levels also in the 
WKB approximation. 



At the singular point t ~ oo of the integrand of Eq. 
||5DJ| Eq. igSJ reduces to 



CONCLUSION 



t dy p 
a dt 



2y p y P -i + e 



(67) 



where c p is a constant satis- 



2c p Cp-i + e. Since at 



whose solution is y p = c p 
fying the algebraic equations c 2 
large p we expect y v -\ — > y p — > y, where y is an exact so- 
lution at t = 0, it means , in view of c p being a constant, 
that should be c p = c p _i = c, that is we are looking for 
a fixed point solution of Eq. (|5T|l which is c p — \fe. As 
in previous paragraphs, the positive sign before the root 
is chosen since the first term in expansion of y p (t) in the 

WKB terms is ik(t). Thus y p {t) ~ ik(t) 

that is y p (oo) ~ +y/e. 

Another singular point of the integrand of Eq. I|50l) 
lays at t ~ 0, where Eq. l(63|l reduces to 



t dy p 
a dt 



vl ~ 2y p y p -i + 



A 



B 

T 



(68) 



We have shown that the quasilinearization method 
(QLM) which approaches solution of the Riccati- 
Schrodinger equation by approximating the nonlinear 
terms by a sequence of the linear ones, and is not based on 
the existence of a smallness parameter, sums the WKB 
series. 

The advantage of the quasilinear approach is that each 
p-th QLM iteration is expressible in a closed integral 
form. We have proved that its expansion in powers of h 
reproduces the structure of the WKB series generating an 
infinite number of the WKB term with 2 P terms of the ex- 
pansion reproduced exactly and a similar number approx- 
imately. As a result one expects that the exact quantiza- 
tion condition (|26Jl with integrand replaced by any QLM 
iterate (|27|) including the first gives more accurate energy 
than the WKB quantization condition l|28|l which is ob- 
tained by substituting into exact quantization condition 
of the WKB expansion up to the first power of h and ne- 
glecting all higher powers of H. We show on many exam- 
ples that it is indeed so and the approximation by QLM 
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iterates leads to exact energies for many well known phys- 
ical potentials with such as the Coulomb, harmonic os- 
cillator, Poschl- Teller, Hulthcn, Hyleraas, Morse, Eckart 
etc. 
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